Quantum correlations as precursors of entanglement 
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We show that for two initially excited qubits, interacting via dipole forces and with a common 
reservoir, entanglement is preceded by the emergence of quantum and classical correlations. After 
a time lag, entanglement finally starts building up, giving rise to a peculiar entangled state, with 
very small classical correlations. Different measures of quantum correlations are discussed, and their 
dynamics are compared and shown to lead to coincident values of these quantifiers for several ranges 
of time. 
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I. INTRODUCTION 

The characterization of entanglement and the elu- 
cidation of its role in quantum computation remain 
formidable challenges, in spite of the conspicuous pres- 
ence of this concept in quantum physics since the fun- 
damental and instigating papers published by Einstein, 
Podolski, and Rosen, as well as Schrodinger, in 1935 
P, 0] ■ Motivation for the understanding of this subtle 
concept is stimulated not only by its fundamental char- 
acter, but also by its applications in quantum informa- 
tion |3| and, much specially, by the perspective that en- 
tanglement could be the key ingredient in the increased 
efficiency of quantum computing compared to classical 
computation for certain quantum algorithms [§-[3 ■ 

In an entangled state, classical and quantum correla- 
tions may coexist. Indeed, for instance, in the singlet 
state |*_) = (| 10) - |01))/V2 there is a perfect classi- 
cal correlation between the first and the second qubit, 
namely if the first is in state |1) the second is in state |0), 
and vice-versa. Several measures of entanglement have 
been proposed, like the concurrence and the negativity 
0- Also, different criteria for the existence of classical 
correlations have been proposed d, Q . In Ref. @ , the 
existence of genuine classical correlations was associated 
to the non- vanishing of n-party correlation functions in- 
volving local observables of the system. Based on this 
definition, it was shown that it is possible to have multi- 
party entangled states with no genuine classical correla- 
tions, as long as the number of parties is larger than two 
1- 

More recently, it has become clear that entanglement 
does not exhaust the realm of quantum correlations. In- 
deed, separable states can exhibit quantum correlations, 
which seem to play a role in the explanation of the power 
of some schemes of quantum computation [Iol - [l2| . Sev- 
eral quantifiers have been proposed for these quantum 
correlations, starting with the work by Ollivier and Zurck 
fl3j , who introduced the "quantum discord" . Intuitively, 
this measure quantifies, in a bipartite system, the mini- 
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mum change in the state of the system and on the infor- 
mation of one of its parts induced by a measurement on 
the other part. For a state with zero quantum discord, it 
is possible to measure any of its parts without changing 
the state of the system. An example would be the state 

pi = ^(|o><o|®|+)(+| + |i)(i|®|-)H)> W 

for which there exist one-dimensional complete projective 
measurements on both the first and second subsystem - 
namely, measurements on the basis {|0), |1)} for the first 
qubit and {|+), |— )} for the second qubit - that do not 
perturb the overall quantum state, nor the state of each 
part. The quantum discord of this state is zero: the state 
has only classical correlations. 
On the other hand, for the state 

P2 = i(|0)(0|®|+)(+| + |l)(l|®|0)(0|) , (2) 

it is clear that any projective measurement performed on 
the second qubit disturbs the state of the first qubit. In 
this case, quantum discord is different from zero. The 
same is true for the maximally-entangled singlet state, 
which is transformed into a statistical mixture upon mea- 
surement of any of the two qubits. 

An interesting implication of this concept is the 
demonstration that vanishing quantum discord is nec- 
essary and sufficient for completely positive maps flij ]. 
that is, it is not necessary to assume, in the usual mas- 
ter equation description, that the initial state of sys- 
tem+environment is a product state. It could be a more 
general separable state, as long as the quantum discord 
vanishes. Quantum discord has been calculated for sev- 
eral families of quantum states and compared with the 
entanglement [15| . Modified versions have been pro- 
posed, with different physical meanings [Lsj . 

Other measures of quantum correlations, involving 
non-orthogonal measurements (positive operator-valued 
measurements - POVM's) and measurements on both 
qubits, have been introduced in the literature [1343 j 
and have been studied and compared with quantum dis- 
cord in different situations: the DQC1 model of mixed- 
state quantum computation [2(|, and within the context 
of accelerated frames [2l|. We present in Sec. HT1 these 
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other quantifiers, as well as a precise definition of quan- 
tum discord. The existence of quantum correlations in 
the absence of entanglement is another subtle trait of 
quantum mechanics, still to be fully understood. 

In this paper, we show that when two qubits inter- 
act with each other through a common environment and 
also through dipole forces, there is a peculiar dynami- 
cal behavior intertwining quantum correlations, classical 
correlations, and entanglement. For two qubits initially 
excited, quantum correlations show up as precursors of 
entanglement, growing up and then shrinking as entan- 
glement belatedly appears. At this moment, classical 
correlations, even though non- vanishing, become prac- 
tically negligible. The dynamics of quantum correlations 
is studied in terms of two distinct quantifiers, the above- 
mentioned quantum discord and the measured-induced 
disturbance (MID), introduced by Luo [l^. In particu- 
lar, we show that the MID coincides with the quantum 
discord for part of the evolution. 

The dynamics of multiparty entangled systems inter- 
acting with independent environments that act on each 
of its parts is quite different from the time-dependent be- 
havior of each individual component of the system. Quite 
generally, the decay of entanglement is non-exponential, 
and it may lead to separability at finite times, before 
each part reaches its final state [22T - I33T ] . The finite- 
time disappearance of entanglement, sometimes called 
"sudden-death of entanglement" [28|, 1331 was experimen- 
tally demonstrated by Almeida et al. [321 ] . The dynamics 
of quantum correlations other then entanglement has also 
been analyzed, both theoretically [U, HH and experimen- 
tally m. 

As the two qubits get closer, at a distance comparable 
to the radiated wavelenght, a different physical situation 
arises. The model of individual and independent envi- 
ronments does not correspond anymore to the physical 
reality: one must consider that the qubits interact with 
the same environment, and furthermore the interaction 
between the qubits, which depends on their physical char- 
acterization, must also be taken into account. This has 
important consequences for the dynamical behavior of 
the system. A common environment may entangle ini- 
tially separable systems, even in the absence of direct 
interaction [37H47| . This can be easily understood from 
the fact that the product state |0) | <8> 1 1) may be expressed 
as a sum of a singlet and a triplet component: 

[01> = |[(|01> - |10» + |01> + |10»] , (3) 

and, since the singlet component does not decay under a 
common environment, a residual entanglement remains 
at asymptotic times [48j |. 

This is not valid anymore for two initially excited 
atoms, since then the singlet component does not ap- 
pear. In this case, one can show that the system remains 
separable throughout its evolution |4§} . The situation 
changes dramatically, however, when dipole interactions 
between the two initially excited atoms arc taken into 



account. In this case, it was shown that entanglement 
does emerge, but only after a time lag pri 149^. 

Here we study the dynamic evolution of classical and 
quantum correlations corresponding to this peculiar be- 
havior of entanglement. As mentioned before, three main 
results emerge from our work: (i) The belated appear- 
ance of entanglement is preceded by the build up of 
strong classical and quantum correlations, which so to 
speak prepare the scenario for the emergence of entan- 
glement; (ii) The interaction between the two initially 
excited qubits, indirectly through the environment and 
directly through dipole forces, leads to an entangled state 
that has extremely small classical correlations; and (iii) 
although stemming from very different definitions, differ- 
ent quantifiers of quantum correlations coincide within 
some ranges of time as the system evolves: this is the 
case for the MID and the quantum discord. 

This paper is organized as follows. In Sec. [TTJ we re- 
view the main quantifiers of quantum and classical cor- 
relations, with special attention to quantum discord and 
the "Measure-Induced Disturbance" (MID) introduced 
by Luo [l8j|. In Sec. IIIII we present the master equation 
that describes the interaction of two qubits with a com- 
mon environment and the main results concerning the 
dynamical behavior of several quantities that have been 
introduced to characterize entanglement and quantum 
correlations: concurrence, quantum mutual information, 
quantum discord, and measure-induced disturbance. We 
also consider in the same section the evolution of the 
classical correlations. The conclusions are presented in 
Sec. IIV1 while the Appendices contain detailed demon- 
strations of some of the assertions made in the article. 



II. QUANTUM AND CLASSICAL 
CORRELATIONS 

A. Quantum discord 

The total classical correlations between two random 
variables X,Y are given by the Classical Mutual Infor- 
mation: 

I(X :Y) = H(X) + H{Y)-H(X,Y), (4) 
where H(X) = - 2jj5psr= x ) logJJp^), is the Shan- 

X 

non entropy, with the probability distributions calcu- 
lated from the joint one, px.Y- Px = ^~]px,Y=y, Py = 

v 

^Px=x,y, and H(X : Y) = ~^2p x . y logp x . y is the 

x x,y 

joint entropy with p xy the probability of both outcomes 
x and y happening. 

Using the Bayes rule, Px\Y=y = Px.y=ylPY=y, where 
Px\Y=y is the conditional probability that the event X 
occurs once the event Y = y has already occurred, the 
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Classical Mutual Information can be expressed equiva- 
lently as: 

J(X :Y) = H(X)-H(X\Y), (5) 

where H(X\Y) = E y PY= y H(X\Y = y) is the con- 
ditional entropy of X given Y, with H(X\Y = y) = 
— ~y^^p x \y logp^iy In this form, it becomes clear that the 

X 

Classical Mutual Information describes the difference in 
the ignorance about the subsystem X before and after 
performing a measurement on subsystem Y . 

Ollivicr and Zurck generalized these two equivalent ex- 
pressions of the Classical Mutual Information to quan- 
tum systems. Equation (QJ is easily generalized by re- 
placing the probability distributions by density matrices 
and the Shannon entropy by the Von Neumann entropy 
S(p) = — Tr(/?log 2 ,o), thus obtaining: 

I(pa,b) = S(p A ) + S(p B ) - S(p A . B ). (6) 

It is not straightforward to generalize Eq. ([5]), since the 
definition of the conditional entropy S(pa\pb) involves 
specifying the state of the subsystem A, knowing the 
state of the subsystem B. In quantum theory, this im- 
plies that a measurement must be applied to subsys- 
tem B. The generalization proposed by Ollivicr and 
Zurek is done by assuming a complete unidimensional 
projective measurement made on system B, correspond- 
ing to the projectors {nf}, such that Hf = 1. 
The state of A after this measurement is implemented 
is PA\nf = n f PA.B^lf /Tr(Hf pa,b), with probability 
Pj = Tr(TY® pa,b)- Then, the conditional entropy is de- 
fined as S {pa\{RB }) = *^^pjS(p A \TiB), and the quantum 
j 

equation corresponding to Eq. (|5|) can be expressed as 

J(PA, B ) { af } = S(pa) - S( PA \{Uf}). (7) 

The quantum discord is the minimum of the difference 
between Eq. ((6]) and Eq. (0: D(pa,b) = mm(J(p^, B ) - 

J{Pa,b){u b })- The quantity J{pa,b){ji b } 1S the infor- 
mation gained about system A when measurements are 
performed on system B. Through the process of min- 
imization over all possible measurements on system B 
we search the measurement that disturbs least the to- 
tal quantum system and allows to obtain as much in- 
formation as possible from system A. Ollivier and Zurek 
demonstrated that 

D{p A , B ) { nf} = O Pa,b = nf PA, B nf , (8) 

3 

that is, a zero quantum discord implies that the com- 
plete quantum system is not disturbed by the measure- 
ment. Furthermore, the information about the system A 
is not perturbed by the measurement of system B. How- 
ever, a nonzero quantum discord implies that the mea- 
surement disturbs the state and part of the information 



about system A that exists in the correlations between 
the subsystems A and B is lost. Equation shows that 
the concept of quantum discord admits a simple physi- 
cal explanation. However, this measure has an annoying 
feature: it is asymmetric under exchange of systems A 
and B, which should not be expected from a quantifier 
of quantum correlations. This is made very clear by go- 
ing back to the example given by Eq. For any com- 
plete one-dimensional projective measurement performed 
on the second qubit, the quantum state of the first one 
is disturbed and therefore information about it is lost. 
The quantum discord for this case is different from zero. 
But if we perform the measurements on the first qubit, 
the quantum state is not perturbed and all the informa- 
tion on the second qubit, initially present in the state, 
is recovered. In this case the quantum discord would 
be zero. This example also highlights the difference be- 
tween separability and classicality, and clearly associates 
the existence of quantum correlations in a separable state 
with the presence of nonorthogonal states in one of the 
subsystems, as already mentioned by [ill ]. As shown in 
fl3j , the quantum discord is always greater or equal than 
zero, and is zero if and only if the state has only classical 
correlations. 



B. Measurement-Induced Disturbance 

Given any bipartite state p, one may associate with 
it, by means of local measurements, another state, inter- 
preted as the classical part of the former. The quantum 
correlations present in the state p are then determined 
by quantifying the difference between these two states. 
Consider any complete set of one-dimensional orthogo- 
nal projections n", acting on each party a and b. 

The state after the measurement is: U(p) = J^(n° (g> 

n^) ( o(H° (g) n^). If for some measurement H(p) = p, then 
the state is called classsical state, otherwise the state is 
truly quantum. Luo demonstrated in fl8j that, if p is 
classical, {n°} , {n^} and {Uf <g> n^} are the eigenpro- 
jectors of p a — Tr^p, p^ = Tv a p and p respectively. This 
implies that the definition of classical states is unam- 
biguous, since there is a unique measurement that leads 

to n(p) = p. 

When a complete set of projective measurements {n^} 
is performed on a system with state described by the 
density matrix p, the entropy of the final state n(p) = 
HipHi is greater or equal than the entropy of the ini- 
tial state, S(p) < S(H(p)), and equality is attained only 
when the projective measurements are the eigenprojec- 
tors of the matrix p [Hoj ] . In order to quantify the quan- 
tum correlations in a quantum state p, Luo [l8j chose 
the measurement n induced by the eigenprojectors of 
the reduced subsystems. With that choice, the reduced 
states remain invariant, and then also the correspond- 
ing entropies: the information on each subsystem is not 
changed and, in that sense, this measurement is the least 
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disturbing. 

The corresponding measure of quantum correlations, 
named "measurement-induced disturbance", is defined 
as: MID( ( o) = I(p)—I(U(p)), where I is the quantum mu- 
tual information and II = II? ® 11$, with II? and 11$ the 
eigenprojectors of the reduced subsystems, p a = T^b(p) 
and pb — Tr Q (p), respectively. In particular when p is a 
pure bipartite state, MID(p) — D(p) = S(p a ): this mea- 
sure coincides with the quantum discord, and is equal to 
the entropy of the reduced system, which is a measure 
of entanglement. The MID presents a great advantage 
with respect to the quantum discord: it is easily calcula- 
ble, since it does not involve any minimization, one has 
only to find the eigenvectors of the density matrices cor- 
responding to the subsystems. However, it is not appli- 
cable in all cases as Wu et al. have remarked: when 
the local density matrices have degenerate eigenvalues, 
the MID is not uniquely defined. 



C. Classical Correlations 

A quantifier for classical correlations was proposed by 
Henderson and Vedral [l7|- The Classical Mutual In- 
formation given by Eq. © was generalized to quantum 
mechanics by replacing the Shannon entropy by the Von 
Neumann entropy and the classical probability distribu- 
tions by density matrices. The corresponding quantifier 
is given by: 



in a bipartite state, they define Q(pab) = I(pab) — 
Ic max {pAB)- When pab is a pure state Q(pab) = 
S(pa) = D(pab) = MID(pab), and in the case of a 
general mixed state, MID(pab) > Q{pab) > D(pab)- 

In this work we refrain from studying the dynamical 
behavior of the quantity Q, because it requires an opti- 
mization over all possible measurements, and an analyt- 
ical expression for it is still unknown. 



D. Entanglement 

For the quantification of entanglement we use the Con- 
currence j52[, defined as: 

C = max{0,A}, (10) 

where: 

A= VAT- V%- V^3- aAL (11) 
Xi being the eigenvalues in decreasing order of the matrix: 

P(<7y ® <Ty)p*(<Ty ® (Ty), (12) 

where p is the two-qubit density matrix, a y is the second 
Pauli matrix and the conjugation is performed in the 
computational basis. Concurrence ranges from 0, which 
corresponds to a separable state, to 1, which corresponds 
to a maximally entangled state. 



Ca(pab) = ma,x(S(pB) 



PiS(p B )), 



(9) 



where p l B = Tr a {AiPABA\)/Tr(AiPABA\) is the state of 
the subsystem B after performing the POVM Ai on A. 

Hamieh et al. showed in [HT| that, for the case of two 
qubits, the POVM that maximizes Eq. ^ is a complete 
set of unidimcnsional orthogonal projectors. With this 
result it is easy to see that I(pab) = D(pab) +C(pab), 
thus the total correlations of the system, quantified by 
the Quantum Mutual Information ©, are separated in 
quantum correlations (measured by the quantum dis- 
cord) and classical correlations. One should note that 
this definition also presents the problem of asymmetry: 
the quantification of classical correlations depends on 
which subsystem is measured. 

Wu et al. [l9| proposed an alternative definition for 
a quantifier of quantum correlations, which is symmetric 
and, as opposed to MID, unique even when the states 
of the subsystems have degenerate eigenvalues. They 
proposed to quantify the classical correlations between 
two systems A and B by performing a POVM locally on 
each subsystem. With the detection records they calcu- 
late then the Classical Mutual Information for the chosen 
POVM. 

They define I Cmax {pAB) as the maximal classical mu- 
tual information available over all choices of possibles 
POVM's. In order to quantify the quantum correlations 



III. DYNAMICS OF CORRELATIONS AND 
ENTANGLEMENT 

A. Theoretical Model 

We consider a system of two identical qubits interact- 
ing with all modes of the electromagnetic field, assumed 
in the vacuum state. The state of each qubit is repre- 
sented in the basis {|e) , \g)} (e = excited, g = ground 
state). The two qubits interact via dipole forces, associ- 
ated to the dipole transition moments p.The total Hamil- 
tonian of the atoms plus the electromagnetic field, in the 
electric dipole approximation, is 



H.c 



where Sf = |e.j) (<?j| and S~ = \gi) (e»| are the lad- 
der operators, Sf = |e») (e« j — \gi) {g%\ is the energy op- 
erator of the ith qubit, u>i are the transition frequen- 
cies (in what follows wc will consider all frequencies 
equal, uji = ojq), at and at the annihilation and cre- 
ation operators corresponding to the field mode ks, with 
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wave vector k, frequency u>k and index of polarization s. 
dks^i) = ( ) %s e r< 1S the coupling constant, 



\2e hV J ks 
?i is the position of the ith qubit, V is the normalization 
volume, and is the unit polarization vector of the 
field. When the length of the system is small compared 
to the radiated wavelength (ujq/c), we may neglect the 
spatial variation of <?£ s (ri), so that, in the rotating- wave 
approximation, the Hamiltonian reduces to [53^ : 



H = huj S z + Y 



u k a\ at 

ks h 



ihYy^-9ks S+ Hs -H.c], 



where S ± and S z are collective spin operators defined by: 

s ± = y s t and sz = Y s ^ 

i i 

This Hamiltonian describes the Dicke model. The 
other limit, when the length of the system is much greater 
than the resonant wavelength, is easily obtained from the 
master equation approach that we considerer in the fol- 
lowing. 

The dynamical evolution of the qubit system is given 
by the following master equation [461 ]: 



d P (r) 

at 



i=l 
2 



if- J 



\ Y r ij (PS?SJ + S+Sjp - 2S- P S+) . 



'>J = 1 



(13) 



where Tu = T are the spontaneous emission rates of the 
qubits, assumed to be identical, and Ti2 = T21 and J7i2 
are respectively the collective damping and the dipole- 
dipole interaction defined by 



r 12 = -Ti 



1 - (A-^ia) 



1-3 (/i.fi 2 ) 



sin(fc ri 2 ) 
korw 

cos(k rij) sin(fc ri2) 



L (V12) 2 



(14) 



and 




cos(fcori 2 ) 

^0^12 

sin(fc ri 2 ) cos(fc ri 2 ) 
{k ri 2 ) 2 (k r 12 ) 3 



(15) 



where fco = ojq/c, t\ 2 = \v\ — r 2 | is the distance between 
the qubits, and p, as before is the dipole transition mo- 
ment. It follows from Eq. (PJ) that 7 = r 12 /r < 1. 

When the distance between the qubits is much greater 
than the resonant wavelenght, it is easy to see that 



ftij — > and Tij 0,i ^ j, so that the master equa- 
tion becomes 



d P (t) 
dt 



-iw ^[5f,p] 

i=l 

1 2 

-TYipStS- + S+S-p - 2S- P S+) . 



which corresponds to the case of independent environ- 
ments. 



B. Results 

In order to study the dynamical evolution of quantum 
and classical correlations, the master equation (|13p was 
solved analytically; we refrain from showing the details of 
the solution here since our results coincide with those in 
Ref. [4|| . From the analytical solution, we calculated the 
relevant quantities for our purposes: concurrence, quan- 
tum discord, MID and classical correlation, for several 
initial states. 

The dynamical evolution of these quantities when both 
qubits are initially excited, interacting with a common 
reservoir and with each other via dipole-dipole interac- 
tion, is shown in Fig. [TJ We see that initially all correla- 
tions, with the exception of the concurrence that remains 
zero, rise to a maximum value and then decrease. One 
should also note that, for a finite period of time, the MID 
and the quantum discord coincide. This can be under- 
stood by using the results in j54| ; for a subset of density 
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FIG. 1: (Color online) Evolution of quantum correlations: The 
thick black line is the quantum mutual information, the thin 
blue line is the quantum discord, the dashed yellow line is the 
measured-induced disturbance, the green dotted line is the classi- 
cal correlation, and the red dottcd-dashed line is the concurrence. 
(r 12 = A/8,7 = 0.8806) 
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matrices with structure X, 



P = 



a lu 

b z 

z b 

uj d 



the quantum discord D{p) is given by: D(p) 
min{Di, -D 2 }, where 



D 1 = 

and 
D 2 



S(/) 
dlog 2 



S{ P AB ) 
d 



b + d 



~ alog 2 
b log, 



a + 6 
6 



61og 2 



b + d 



a + b 
(16) 



S^) - S( P AB ) - |(1 + a) log 2 (1 + a)) 

_!(l-a)]og a (i(l-a)), (17) 

with a 2 = (a - rf) 2 + 4\z + u\ 2 . 

The density matrix describing the dynamical evolu- 
tion of the two-qubit system here considered, when both 
qubits arc initially excited, has the same form as the one 
in Eq. (TITjI) . Indeed, we get 



P(r) 



where 



/ a(r) 

6(r) c(r) 

c(r) 6(r) 

\ l-a(r)-2&(r) 



(18) 



a(r) 
6(r) 

c(r) 



with 



s~ 2t 

he- 2r + e-( 1 -^) r ](l- 7 ) 

2(1 + 7) 
[-e- 2r + e-( 1+ T) r ](l + 7 ) 

2(1-7) 
[-e-^ + e-^-^Kl-T) 
2(1 + 7) 
[-e- 2r + e-( 1+ T) r ](l + 7 ) 

2(1 = 7) ' 



r = rt, 7 = r 12 /r. 



(19) 



(20) 



As a function of the renormalized time r, these expres- 
sions depend on a single parameter 7. It is easy to show 
from thie above equations that the analytic expression 
for the MID coincides with the expression for D\ (see 
Appendix A), so whenever D\ < £> 2 , the MID coincides 
with the quantum discord. This condition is fulfilled for 
a finite time domain, that depends on the parameter 7 
as shown in Fig. 2J 



At a later time t e , the system becomes entangled and 
the concurrence starts increasing. This feature was previ- 
ously found by Ficek and Tanas in J4(| H^] . They showed 
that, for two qubits initially in their excited states, under 
the conditions considered in this paper (interacting with 
a common environment and through dipole forces), and 
for small distance ri 2 between them (as compared to the 
resonance wavelength), entanglement appears only after 
some finite period of time. This feature is independent 
of the values of the parameters, as shown in Appendix B. 
After the time interval where MID and quantum discord 
match, the numerical value of the quantum discord re- 
mains close to the MID, but always a little lower. On the 
other hand, after the emergence of entanglement, the dif- 
ference between MID and Concurrence continuously de- 
creases, going to zero when both vanish at infinite times. 

Classical correlations are smaller than quantum corre- 
lations for most of the dynamic evolution of the state, 
a result that was previously found in (55j . After they 
reach their maximum value, they go to zero asymptoti- 
cally. However, for most of the time of the evolution, they 
have a very small value, about a hundred times smaller 
than the concurrence. This behavior is quite general: we 
have demonstrated it for several values of the parameter 
Ti 2 , as shown in Appendix C. In Ref. @, it was shown 
that only with three or more qubits can one have states 
with genuine quantum correlations but no classical cor- 
relations. Our example does not contradict that general 
result, since for the two-qubit system here considered the 
classical correlations, as defined in Q, are indeed differ- 
ent from zero, but it also shows that the same correla- 
tions, even though not vanishing, may be much smaller 
than the quantum correlations, even for a two-qubit sys- 
tem. 

Ficek and Tanas also showed in Ref. [49| that, for 
the Dicke model (evolution of two two-level atoms in a 
common reservoir with no dipolc-dipolc interaction) , two 
qubits initially excited do not get entangled. A question 
naturally arises as to whether that is the case for the 
other quantum correlations in the Dicke model. This 
question is answered by taking the limit ri 2 — > in the 
expressions for Ti 2 and f2i 2 . In that case, Ti 2 — > T, so 
that 7 -> 1, and fii 2 -> [3r/4(fc r 12 ) 3 ] [1 - 3(p.ri 2 ) 2 }. Al- 
though rii2 diverges in this limit, Eq. (jT9")l shows that the 
evolution of the initial quantum state considered here is 
independent of the dipolc-dipolc interaction. Therefore, 
in this limit we do recover, for the considered initial state, 
the evolution predicted by the Dicke model. The density 
matrix for this model is easily obtained from Eqs. (fTJ 
and (TlT)|) taking the limit 7 — > 1: 



P(r) 



where 



/ o(r) 

6(r) 6(t) 

6(r) 6(t) 

V 



\ 



(21) 



1 - a(r) - 26(t) / 



a(r) 
b(r) 



er 2T T. 



(22) 
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FIG. 2: (Color online) Evolution of the quantum correlations for 
the Dicke model: the thick black line is the Quantum Mutual Infor- 
mation, the thin blue line is the quantum discord, the dashed yellow 
line is the MID, the green dotted line is the Classical Correlation 
and the red dotted-dashed line is the Concurrence. 



Figure [2] displays the evolution of all the above- 
mentioned correlations, within the framework of the 
Dicke model, when the two qubits are initially excited. 
In this case, the initial behavior of classical and quan- 
tum correlations (except entanglement) is the same as 
in the model that includes dipole-dipolc interaction: an 
initial period of growth to a maximum and then a decay. 
However, one does not have, as in the former situation, 
a revival of these correlations. 

The lack of entanglement for an initially doubly ex- 
cited two-qubit state |ee) in the absence of dipole-dipolc 
interaction stems from the fact that the singlet state 
^7f (I e d) ~ l<? e )) does not show up in the evolution, as dis- 
cussed by Ficek and Tanas [46| : this is a decoherence-free 
state in the Dicke model, and its constant contribution 
leads to asymptotic entanglement when the initial state 
is | eg), since this state can be written as a sum of the 
anti-symmetric singlet ^(l e 5) — \g e )) an( i the symmet- 
ric state 4^ (| eg) + \ge)). The dipole-dipole interaction, 
on the other hand, leads to the emergence of the singlet 
state for the initial state |ee), to the survival of quantum 
correlations, and to the belated appearance of entangle- 
ment. Indeed, the analytical expressions for the MID and 
quantum discord show that, for r >> 1, the decay time 
of the symmetric state, both MID and quantum discord 
decay as e -( - 1-7 - |T , which is the decay law for the pop- 
ulation of the singlet state. This is shown in Figure G2 
which makes clear the essential role of the singlet state in 
explaining the remarkable behavior of the entanglement 
and the quantum correlations in the presence of dipole- 
dipole interactions. 

When the distance between the qubits r 12 becomes 
much larger than the resonant wavelength, we recover the 
expected behavior for two initially excited qubits evolv- 
ing in independent environments: the state remains sepa- 
rable and no correlations are created between the qubits. 
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FIG. 3: (Color online) Evolution of the quantum correlations 
when the dipole-dipolc interaction is included: The green solid line 
is the quantum discord, the blue dashed line is the MID, the red 
thick dashed-dotted line is the Concurrence, the brown dotted line 
is the population of the symmetric state and the orange thick line 
is the population of the antisymmetric state. The inset shows the 
zone inside the ellipse. (r\2 = A/4,7 = 0.9460) 



IV. CONCLUSIONS 

In this work, we have studied the subtle dynamics of 
quantum and classical correlations for two qubits, ini- 
tially in a pure product state, coupled through dipolc 
forces and interacting with all modes of the electromag- 
netic reservoir. When the two qubits are initially ex- 
cited, entanglement has a peculiar behavior [H, it 
remains zero for a finite-time interval, and then it builds 
up. Our work is aimed at clarifying what happens dur- 
ing the "dormant" time, that is, what kind of dynamic 
changes prepare the system for the late appearance of 
entanglement. With this aim, we have analyzed the dy- 
namic behavior of quantifiers of quantum and classical 
correlations. We show that, for two excited qubits, the 
reservoir creates initially classical correlations and quan- 
tum correlations between the qubits, which remain in 
a separable state. These correlations evolve from zero 
to a maximum value, and then decay. This overall be- 
havior does not depend on the presence of dipole-dipole 
interactions. However, this interaction is fundamental 
for keeping up the quantum correlations and building up 
entanglement after a time lag. It delays the decay of 
quantum correlations, thus allowing for the build up of 
entanglement. 

Therefore, in the presence of dipole-dipole interactions, 
quantum correlations can be considered as precursors of 
entanglement. In the absence of these interactions, they 
build up and decay to zero, the state remaining separable 
for all times. The dipole-dipole interaction helps preserv- 
ing quantum correlations, and this seems to fire up en- 
tanglement in this case. While entanglement is still zero, 
there is a time span for which MID and quantum discord 
coincide. Even before entanglement appears, classical 
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correlations become negligible, and remain so throughout 
the evolution. Therefore, the very evolution of the sys- 
tem generates a peculiar entangled state, with very small 
classical correlations. Even though there is still consid- 
erable controversy over the proper definition of quantum 
correlations and their role in quantum computation, their 
dynamics under the action of the environment seems to 
be, in the present context, intimately related to the gen- 
eration of entanglement. Further studies in this direction 
might help to elucidate this subtle dynamic connection 
between quantum correlations and entanglement. 

V. APPENDIX 




A. MID equals quantum discord 

The density matrix corresponding to the initial state is 
given by Eq. (|T8"j) . To calculate the MID, we apply projec- 
tive measurements on the density matrix, corresponding 
to the eigenvectors of the reduced density matrices. Af- 
ter these measurements, which amount to keeping only 
the diagonal elements in Eq. (|18[) . the density matrix is: 

/ a(r) \ 

r \ _ o b(r) 

PmeasKJ) — Q &(t) Q 

V 1 - a(r) - 26(t) / 

Then 

MID(p) = -S(p) + S(p meas ) = -26(r)log 2 [6(r)] 
+ [6(r)-c(r)]Iog a [6(r)-c(r)] 
+ [b(r)+c(T)}log 2 [b(T)+c(r)}, (23) 

expression that coincides with Dt, Eq. (flH]). through- 
out the evolution. According to |54j the quantum dis- 
cord is the minimum of {Di,!^}- Therefore, whenever 
Di < Z?2, the MID coincides with the quantum dis- 
cord. The expressions for D\ and D 2 do not allow get- 
ting an analytical solution for the interval of time where 
D\ < D 2 - Fig. [4] exhibits the results of numerical cal- 
culations made for several values of 7, displaying time 
intervals (reduced to a single point for small 7) where 
D 1 < D 2 . 

As mentioned above, in the calculation of the MID 
there is an ambiguity when the reduced matrices have 
degenerate eigenvalues, or equivalently, are multiples of 
the identity matrix. In this case, any two orthogonal 
vectors are eigenvectors of the reduced matrices, so there 
are infinite ways to choose the local projective measure- 
ments, which give rise to different values of the MID. 
Here the reduced density matrices are: 

, , ( . ( o(t) + b{r) \ 

P«lr)=P»<r)=\ l-a(r)-6(r)J ' 

and the mathematical condition for the reduced matrices 
to become proportional to the identity is: 

o(r) + 6(r) = i (24) 



FIG. 4: (Color online) Evolution of the expressions D\ (blue line) 
and D2 (red dashed line) for different values of the parameter 7. 



with 

o(r) + fe(r) = e~ 2T + £ {[cosher) - e~ r ] 
(I-7 2 ) 

x (1 + 7 2 ) — 7 sinh(7T)} . 

This expression decreases monotonically from its initial 
unitary value to zero, as r — > 00 - see Fig. [5] There- 
fore, there can only be a single instant of time for which 
Eq. (f2~4")l holds. At this point, the MID is not uniquely 
defined and depends on the eigenvectors chosen to make 
the measurements. However, since, for all other times, 
the eigenvectors are uniquely defined and coincide with 
the computational basis, it is natural to choose the same 
basis for the degeneracy points, in which case the MID 
does not present discontinuities. 




12 3 4 

T 



FIG. 5: (Color online) Plots of a(t) + b(t) for 7 = O.f (blue 
dashed line), 0.4 (red full line), 0.7 (green dashed-dotted line) and 
0.9 (black dotted line). 
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B. Concurrence 

The concurrence C is expressed in terms of the square 
roots of the eigenvalues of the matrix in Eq. ([12")) : 



Wa(r) [l-a(r)-26(r)], y/a(r) [1 - a(r) - 2b(r)}, 

|6(r)-c(r)|,|6(r) + c(r)|}. (25) 

First we will show that b(r) > c(r). This is immediate 
since, from Eq. ([T9|) . 



b(r) - c(r) 



1+7 



(l- 7 )(e 7T 



We distinguish three different cases: (1) 
\J 'a(r) [1 — a(r) — 26(r)] is the largest value in 
Eq. (f2g]l. in which case C = max{0, -26(r)}, (2) 
|6(t)+c(t)| is the largest value, in which case 



C = 



max 



0, 2 (c(t) - y fl (T) [1 - a(r) - 26(r)]) 
and (3) |6(r) — c(r)| is the largest value, in which case 



C = 



0, -2 (c(t) + y/a(r) [1 - a(r) - 2b(r)]) 



In the first case, it follows from the expression for 6(r) 
in Eq. (|X9|) and from the inequality cosher) — e~ r > 
sinh(7r) that 6(t) > for all r. Therefore, the concur- 
rence is zero in this case. 

In the second case, we must investigate the behav- 
ior of c(r) — ^/a(r)[l — a(r) — 2b(r)]. First, we note 
that this case is realized only in the interval where 
c(t) > 0, so that c(r) 2 < 6(r)c(r). We show now 
that 6(t)c(t) < a(Y)(l — a(r) — 26(r)), which implies 
that c(t) < ^/a(r)[l — o(t) — 26(r)]. This is equivalent 



to proving that /(r) = a(r) + 2b(r) + 



Ht)c(t) 
o(r) 



< 1. 



Since /(0) = 1, the equality is verified at r = 0. Fur- 
thermore, it is easy to show that /'(r) < for all r, 
independently of the value of the parameter. There- 
fore, /(t) is a monotonically decreasing function of r 
with /(0) = 1, so /(t) < 1 for all t, implying that 
c(t) < y/a(r)[l — a(r) — 2b(r)]. Thus, we conclude that 
for this second case the concurrence is also zero. 

In the third case, we must investigate the behavior 
of g(r) = c(t) + y/a(r)[l — a(r) — 2b(r)]. Because the 
equation g(r) = contains the time in a non-algebraic 
expression, it is not possible to determine analytically 
the time when the state becomes entangled; however we 
can show, through some simple mathematical arguments, 
that, whatever the value of the parameter, this time al- 
ways exists. This requires to show that there is a time in- 
terval for which two conditions are satisfied: (i) g(r) < 0; 
and (ii) \b(r) — c(r)| is the largest value in Eq. (|25|) 
Condition (i) is demonstrated by noting that 



lim o(t) = lim 



lim - 




= + , and 



lim g{r) 



Then between zero 



and infinite there exists an instant of time r e for which 
g(r) = 0. This implies the existence of a time interval 
from r e to infinity where g(r) is negative, which would 
imply entanglement, as long as condition (ii) is satisfied 
within at least part of this time interval. 

We show now that condition (ii) holds for suffi- 
ciently large r. This follows from lim c(t) = 

" ^ --^-t) 7 - = 0", which implies that 



lim 



>+oo 2(1+7) 

c(r) < for sufficiently large times. Therefore, for times 
sufficiently large, one has g(r) < and \b(r) — c(r)| > 
\b(r) + c(t)|. It remains to show that, for sufficiently 
large times, \b(r) — c(r)\ > \/a(r) [1 - a(r) - 26(t)]. 
This is shown by noting that there is a time t' c such that, 

(1+7) 



for all t>t' s , e 7T 



> 



-, which is immediate. 



(1-7)' 

Therefore, for sufficiently large times, both conditions (i) 
and (ii) above are fulfilled, and the state becomes entan- 
gled. 



C. Classical Correlations 

As mentioned above, classical correlations are much 
smaller than the concurrence for most of the evolution. 
Depending on the analytical expression of the quantum 
discord {D\ or D 2 ), we have two expressions for the clas- 
sical correlations, each valid in a different time interval: 

CCi(r) = - 2[o(r) + b(r)} log 2 [a(r) + 6(r)]- 

2[l-o(r)-6(r)]log 3 [l-a(r)-6(r)]+ 
[1 - o(t) - 26(r)] log 2 [l - a{r) - 26(r)]+ 
a(r) log 2 a(r) + 26(r) log 2 6(r) 



(26) 



and, 
CC 2 {t) 



= -[o(r)+6(r)]log 3 [a(r) + 6(r)]+ 
a(r) log 2 a(r) + [6(t) - c(r)] log 2 [6(r) - c(r)]- 
[6(r) + c(r)] log 2 [6(r) + c(r)]+ 
[1 - a(r) - 26(r)] Iog 2 [l - a(r) - 26(r)]+ 
[o(r) + 6(r)]log 2 [l-a(r)-6(r)]+ 



a) log 2 



a) 



(1-a) 



(27) 



where a(r), 6(r), c(r) are given in Eq. (|19p and a = 
[2(a(r) + 6(r))-l] 2 +4|c(r)| 2 . Figure © displays the clas- 
sical correlations and the concurrence for several values 
of the parameter 7. Typically, classical correlations are 
of the order of 1% of the concurrence. 
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FIG. 6: (Color online) Evolution of concurrence (red dashed line) 
and classical correlations (green line) for several values of the pa- 
rameter 7. 
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